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Electromagnetic Casimir stresses are of relevance to many technologies based on mesoscopic de- 
vices such as MEMS embedded in dielectric media, Casimir induced friction in nano-machinery, 
micro-fluidics and molecular electronics. Computation of such stresses based on cavity QED gen- 
erally require numerical analysis based on a regularization process. A new scheme is described 
that has the potential for wide applicability to systems involving realistic inhomogeneous media. 
From a knowledge of the spectrum of the stationary modes of the electromagnetic field the scheme 
is illustrated by estimating numerically the Casimir stress on opposite faces of a pair of perfectly 
conducting planes separated by a vacuum and the change in this result when the region between 
the plates is filled with an incompressible inhomogeneous non-dispersive dielectric. 
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INTRODUCTION 



Electromagnetic interactions between electrically neu- 
tral isolated polarizable sources such as atoms or 
molecules are often referred to as Van der Waals forces. 
In a quantum field theoretical description they give rise 
to Casimir forces, particularly when some of the sources 
are replaced by a neutral continuum. Such a continuum 
may be restricted to interfaces between different regions 
of space and some regions may contain polarizable media 
with conducting or dispersive properties. In such cases 
one is confronted with the problem of calculating quan- 
tum induced stresses in such regions and the resulting 
pressures on the surrounding interfaces [H, Q. In sit- 
uations where the sources are idealized to occupy per- 
fectly conducting surfaces, it is possible to estimate the 
induced Casimir stresses on the surfaces in terms of quan- 
tum fluctuations of the electromagnetic field in the vac- 
uum . For a pair of approximately parallel conducting 
planes such integrated stresses have been detected ex- 
perimentally and their dependence on the separation be- 
tween the planes measured. However, the analytic deriva- 
tion of vacuum induced Casimir stresses on conducting 
surfaces with non-planar geometries is notoriously dif- 
ficult [3| to ascertain with confidence. These difficul- 
ties are compounded in situations where such surfaces 
bound dielectric media that may be dispersive d, 0] or 
contain inhomogeneous 0, Q magnetic and electric sus- 
ceptibilities Furthermore, the conceptual basis 
on which such calculations are expected to be reliable 
depends on whether it is reasonable to treat polarizable 
systems as a continuum when quantum effects become 
significant. However, one expects a continuum model of 
rigid dielectric media to be reasonable in mesoscopic sys- 
tems where bounding geometries cannot resolve molecu- 
lar detail. In such cases, the response of the medium to 
electromagnetic fluctuations is given in terms of piece- 



wise smooth susceptibility tensors with components that 
may depend non-linearly and non-locally on space and 
time. For media with linear piecewise homogeneous non- 
conducting response functions, Lifshitz |lll4l3t developed 
a phenomenological scheme based on the analytic prop- 
erties of Green tensors. Since its inception, this the- 
ory has not been signiflcantly reflned to deal with more 
general physical systems that are now of relevance in a 
number of modern technologies. These include the influ- 
ence of Casimir stresses due to complicated geometries 
in MEMS devices embedded in inhomogeneous or non- 
linear dielectric media, Casimir induced friction [3| in 
nano-machinery, micro-fluidic and molecular electronic 
devices. Perhaps the most significant feature of the Lif- 
shitz theory limiting its applicability to such systems is 
its reliance on a detailed knowledge of a Green tensor 
(and its analytic structure) leading to a viable regular- 
ization scheme. All attempts to apply the methods of 
cavity QED to mesoscopic systems containing dielectrics 
with inhomogeneous permittivities also rely on a knowl- 
edge of the quantum Hamiltonian of the electromagnetic 
field in the medium in order to calculate finite quantum 
expectation values of "observables" that can be compared 
with experiment. Since any quantum field is an infinite 
dimensional dynamical system, such values need to be de- 
termined by a regularization process that discards unob- 
servable self- forces between sources 15| . When canonical 



dynamical variables can be chosen so that the Hamilto- 
nian for the electromagnetic field has a discrete angular 
frequency spectrum {a;^} and the same structure as the 
Hamiltonian describing an infinite number of simple har- 
monic oscillators at each point in the medium, the regu- 
larization of the (zero temperature) ground state energy 
ift^^cj^ is often defined by continuing to s = — 1 the 
function ^h({s) where 
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with s > So for some Sq that renders the integral con- 
vergent. However, there are few cases in which the sum 



chosen so that each integral 



can be performed analytically and so the an- 



alytic continuation becomes difficult and recourse to nu- 
merics is often inevitable [l6| . Alternative regulariza- 
tion schemes bypass the ground state energy and 
employ point-splitting techniques on components of the 
electromagnetic stress-tensor prior to numerical analy- 
sis. Such approaches are not universally applicable [1| 
and where used often employ the truncation of an infi- 
nite series whose radius and rate of convergence is rarely 
known. 



LAURENT REGULARIZATION 

It is the aim of this letter to offer a more robust nu- 
merical algorithm that has the potential for wider appli- 
cability than existing regularization schemes. It will be 
illustrated by estimating numerically the original Casimir 
stress on opposite faces of a pair of perfectly conducting 
planes separated by a vacuum and the change in this re- 
sult when the region between the plates is filled with an 
inhomogeneous non-dispersive dielectric. It requires for 
its implementation a knowledge of the spectrum of the 
stationary modes of the electromagnetic field between the 
plates. Such a spectrum will in general arise as the infi- 
nite number of roots of a system of transcendental equa- 
tions obtained from imposing the appropriate boundary 
conditions for the modes of the electromagnetic field. 
This involves solving a classical boundary value problem 
for stationary states. For a simply connected bounded 
rigid rectangular cavity with perfectly conducting walls 
such modes can be classified as either TE or TM with, in 
general, different spectra. They are used to construct a 
Fock space of electromagnetic modes in a gauge in which 
the quantum Hamiltonian is quadratic in creation and 
annihilation operators for each mode. For a stationary 
system there is no ambiguity in the classical Maxwell 
stress tensor for the electromagnetic field in a dielectric, 
so one can write down the ground state expectation value 
of any component of this tensor in the medium. Evalu- 
ated on one side of any plane surface this yields, after 
regularization, the quantum stress on that side. In an 
inhomogeneous dielectric, such stresses will in general be 
different on different faces of the cavity. The simplest re- 
alistic situation is to suppose that the permittivity varies 
smoothly in only one direction in a rigid 3-dimensional 
rectangular box with perfectly conducting boundaries. If 
the dimensions of the box are such that two opposite 
pairs of end planes have much smaller areas than the 
pair separated by the direction of the inhomogeneity it 
can be shown that the difference in end pressures in this 
direction requires the regularization of a couple of dou- 
ble integrals, each of the form y) dx dy. Each 
J^(a;, y) is determined by a spectrum generating function 
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T{x, y) dx dy 



regarded as a function in the complex s-plane is analytic 
in an annular region centered on s = 0. As such, it admits 
a representation as a Laurent expansion in this domain. 
The regularized value of J{Q) is then defined to be the 
term cg in this Laurent expansion that is independent 
of s and corresponds to discarding the principal part of 
the Laurent scries before taking the limit as s tends to 
zero. The problem is how to determine numerically cq 
from a numerical computation oi J^{s) when neither the 
principle part of its Laurent expansion about s = nor 
its domain of convergence is known a priori. When J^{s) 
is meromorphic with a pole of order A^i at s = 
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for real s G [e^ , s/j] , finite negative integer iVi and con- 
stants < Eg 1, Sfl > 0. The algorithm for estimating 
Co proceeds first by discretizing the range s e [eg , sn] for 
some sn to generate the set 5 = {sj | 1 < j < J} and 
then evaluating J7(s) numerically at S to generate the set 
X = {Jj = J{sj) \1 < j < J}. Let A denote a matrix 
where each element represents a truncated Laurent series 
of the form: 
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ni < 712 (2) 



for some positive integer 772 and negative integer 771. To 
effect a numerical fit of J'(s) to the "appropriate" Lau- 
rent series, one first determines the matrix elements in 
any sub-matrix of A by fitting the data X, S to each 
(s) by linear regression. Thus, for all integer ranges 
from 77i to 772 with A'^i < ni < — 1 and 1 < 772 < iV2 one 
may calculate the set {c„ (771, 772)} associated with each 
matrix element. We seek criteria such that cq in ([T]) is 
approximated by some 00(771,772) in A, for an "optimal" 
choice of the integers A^i , A^2 defining the size of the sub- 
matrix. The strategy is then to prune the principal part 
of each matrix element by filtering out of each truncated 
Laurent series those terms with coefficients €,1(771,772) 
with 77 < that satisfy 
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Mn{ni,n2) 



for some tolerance < ec <C 1 where the average 
Mn{ni,n2) = j^Y^^Zni c„(j, 772). After this filtering 
process, one has a sequence of truncated Laurent series 
{L"j(s)} whose principal parts contain only the coeffi- 
cients 
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Let Cj\f{ni,n2) denote the coefficient of the most singular 
term in each {i5jj(s)} where by hypothesis A/" < 0. By 
examination of this coefficient for ah -/Vi < ni < — 1 and 
1 < "■2 ^ -^2 in the selected sub-matrix, one may discover 
a new sub-matrix whose elements contain truncated Lau- 
rent series with the same fixed value of M. Next, for each 
positive integer n2, define the data set 

iin2)^{Jjin2)\l<j <J} 

where J} (712) = J} — Cj^{N ,n2)s^ and use linear re- 
gression again to fit this to a sequence of new truncated 
Laurent scries 
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with 1 <n2,n2 < N2. If, for such a fixed J\f, one connects 
the points obtained by plotting Co (A/", ^2 , "■2 ) against n2 
for each 712 in the range 1 < ?i2 < N2, curves are gen- 
erated that attain turning points (or asymptotes) with 
ordinates at co(A/', 712,712) in close proximity. As a result, 
an estimate of the required Casimir coefficient co(A/') is 
defined to be J2no=i co(A^, "2, '^2)- 



APPLICATIONS 

If -F'(z) is meromorphic, f{z) entire and well behaved 
in the complex z-plane where < arg(z) < 7r/2 and {zr} 
denotes simple roots of the equation F{z) = 0, then by 
Cauchy's theorem 
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where 
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For the empty rectangular box with sides of length 
Lx, Ly, Lz, the TE and TM spectra in the box are de- 
generate with: 
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spectrum generating function is 



F{uj, k) = sin L 
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Then, for L^^Ly (i.e. parallel plates) and f{z) = 

Hz/2, the TE Casimir energy {£'^^) is determined from a 
Laurent expansion of 



£;™(s) 



where 



X{s) 



coth(r) dr. 
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For s > 0, the integral ^ can be performed analytically 
and for small s determines all coefficients in the Laurent 
series 



I(s) = ^ C„i 



The TE Casimir energy is then taken to be 

LxLyhc 
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in this case, the value of co can be determined analyt- 
ically since the integral can be written in terms of the 
third derivative of the polygamma function 4'(3,x) as 

The Laurent expansion of this expression gives the exact 
value Co = 7r'*/360. If the integral ^ is evaluated numer- 
ically to 7 significant figures for a range of s in the vicinity 
of the origin, one obtains from our algorithm good agree- 
ment with this analytic result. It agrees with the value 
determined by Riemann ^ function rcgularization within 
1.2%. 
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FIG. 1. For data based on ([3|, the algorithm determines 

Af — —4 and a collection of [x,y] points [n2, Co(A/', n2, 712)] 
with 1 < n2 < 8. The 8 curves obtained by joining these 
points are indistinguishable in this figure and the average of 
co(A/', n2, 712) at their turning points yields a Casimir coeffi- 
cient co{Af) ~ 0.27281. The horizontal dotted line indicates 
the value of the Casimir coefficient that determines the ob- 
served Casimir attractive pressure between perfectly conduct- 
ing plates separated by the vacuum. 

The total (TE+TM) Casimir force per unit area on each 
face is attractive with magnitude: 
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and hence the force difference between the plates is zero. 
Suppose now that the box is filled with an incompressible 
dielectric with inhomogeneous permittivity 

e{x,y,z) = cyiY>{az/L^), < z < 

for some real inhomogeneity parameter a. A quantiza- 
tion of the electromagnetic field in the box [10| can be 
performed in a gauge where the vector potential A sat- 
isfies the condition V • (eA) = 0. For r e {TE,TM}, the 
regularized force differences (AJ^^'^^ (cr)) , derived from the 
quantum expectation value of relevant components of the 
electromagnetic stress-energy-momentum tensor between 
the faces at z = and z = Lz contributed by the TE and 
TM modes in the dielectric is, for L^, Ly 3> L^, derived 
from a Laurent s-cxpansion of 

_ AT^-\s,a) ^ f°°^^-s^a^ r A^^\iy,v,o)e-'y dy. (4) 
•^0 Jo Jo 

Here cr = e"/^ is a positive inhomogeneity parameter, 
To = hca^LxLy/GAir'^L'^ and 

A^^) (z, z^, a) ^ In (i^*^' (z, cr)) 

i^™(z,i/,a) = J^{z)Y^{az) ~ J^{az)Y^{z) 

F™iz, V, a) = J^{z) %{az) - J^{az) %{z), 

where ^ = \J -\-\ and, for any Bessel function Qy^ 
Qy{z) = zQl{z) + Qi/(z). The above double integrals in- 
volve the integration of products of modified Bessel func- 
tions with respect to both order and argument and have 
resisted analytical evaluation. For s > they can how- 
ever be calculated numerically (on a laptop) for various 
cr and fitted by regression to a truncated Laurent ex- 
pansion in the vicinity of s = 0. Our algorithm clearly 
establishes that the principal part of the truncated Lau- 
rent series corresponds to a pole of order 4 (as is the 
case for a = when e is the permittivity of the vac- 
uum). It also determines the Casimir coefficients and 
hence the regularized stress differences between the plate 
faces for both the TE and TM modes. For example, 
with cr = 8/27 the TM modes contribute a force differ- 
ence 3.54704 X 10-28 L^Ly/L\ Newtons. This is 0.27282 
times the value of the total Casimir force on either plate 
in the vacuum. Similarly the TE modes contribute a 
force difference 3.46159 x lO'^s L.^Ly/Li Newtons. This 
is 0.26625 times the value of the total Casimir force on 
either plate in the vacuum. The similarity of these mode 
contributions to each other is somewhat surprising given 
the differences in the structure of the TE and TM spec- 
trum generators. 

SUMMARY 

A robust numerical scheme for regularizing the quan- 
tum electromagnetic stresses in an inhomogeneous di- 
electric between conducting plates has been described. 
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Positive degree of truncated La\irent aeries, 712 

FIG. 2. For data based on Q with o = 8/27 the algo- 
rithm determines M = —4 and a collection of [x, y] points 
[n2,Cj''(A/',n2,n2)] and [n2, (^"(A/", n2, n2)] with 1 < n2 < 8. 
The TE curves obtained by joining the TE points are indis- 
tinguishable in this figure and the average of Co^(A/', fi2, 712) 
at their turning points yields a Casimir coefficient c™(A/') = 
0.19744. The same is true for the TM points with c™(A/') = 
0.20231. Each such Casimir coefficient contributes to the 
total regularized force difference J-q [cq^{N) -|-c™(A/')), be- 
tween opposite z-faces of a pair of perfectly conducting plates 
separated by an inhomogeneous dielectric with permittivity 
e{x,y,z) = eo exp{az/Lz) where a = exp(a/2). 



We believe that it has much wider applicability to more 
general systems such as those mentioned in the intro- 
duction. Although the determination of electromagnetic 
cavity modes in more complex geometries and media is 
in general non-trivial, once this hurdle is overcome the 
quantization program can proceed (with possible ther- 
mal corrections [l^, [2^ ) and the regularization method 
outlined here is then straightforward. However, any reg- 
ularization scheme involving fields in media ultimately 
depends on the viability of the methods of field quanti- 
zation in a non-dynamic background. It is therefore of 
paramount importance to verify precise results of such a 
scheme by experiment. 

Any confined inhomogeneous material dielectric will 
sustain stresses induced by electromagnetic quantum 
fluctuations if the confining domain is rigid. If the 
medium remains static, such stresses induce mechanical 
(elastic) stresses in the dielectric to maintain equilibrium. 
Unlike similarly induced classical stresses by the classical 
gravitational field in the laboratory (that vary with the 
orientation of the dielectric) , the quantum induced elec- 
tromagnetic stresses are permanent. In principle, they 
could be detected experimentally by noting the variation 
of the induced stress field within the dielectric with vari- 
ations of the permittivity inhomogeneities. Such vari- 
ations might be detected using photo-elastic effects on 
the polarization of light passing through a transparent 
medium. 
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